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Thermal Memory in Self-Assembled Collagen Fibril Networks 
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ABSTRACT Collagen fibrils form extracellular networks that regulate cell functions and provide mechanical strength to tissues. 
Collagen fibrillogenesis is an entropy-driven process promoted by warming and reversed by cooling. Here, we investigate the 
influence of noncovalent interactions mediated by the collagen triple helix on fibril stability. We measure the kinetics of cold- 
induced disassembly of fibrils formed from purified collagen I using turbimetry, probe the fibril morphology by atomic force 
microscopy, and measure the network connectivity by confocal microscopy and rheometry. We demonstrate that collagen fibrils 
disassemble by subunit release from their sides as well as their ends, with complex kinetics involving an initial fast release fol¬ 
lowed by a slow release. Surprisingly, the fibrils are gradually stabilized over time, leading to thermal memory. This dynamic 
stabilization may reflect structural plasticity of the collagen fibrils arising from their complex structure. In addition, we propose 
that the polymeric nature of collagen monomers may lead to slow kinetics of subunit desorption from the fibril surface. Dynamic 
stabilization of fibrils may be relevant in the initial stages of collagen assembly during embryogenesis, fibrosis, and wound heal¬ 
ing. Moreover, our results are relevant for tissue repair and drug delivery applications, where it is crucial to control fibril stability. 


INTRODUCTION 

Collagens are the most abundant structural proteins in the 
extracellular matrix of vertebrates. They form filamentous 
frameworks that provide mechanical strength to connective 
tissues such as tendon, skin, and bone (1). Moreover, 
collagen plays a key role in the regulation of important 
cell functions such as migration and differentiation (2). 
The main fibril-forming collagen in mammals is type I 
collagen, which consists of two identical a 1(1) chains and 
a ft'2(l) chain intertwined into a rod-like triple helix. In vivo, 
collagen self-assembles into fibrils upon enzymatic conver¬ 
sion of soluble procollagen precursors into insoluble tropo- 
collagen monomers (3). Tropocollagen has a triple helical 
region of -1000 amino acid residues flanked on both sides 
by short, nonhelical regions known as the N- and C-telopep- 
tides (Fig. 1 A, see figure caption for definition of length 
scales). The monomers are very long (300 nm) and thin 
(1.5 nm) and behave in solution as flexible polymers with 
a persistence length of -10 nm (4,5) (Fig. 1 B). They self- 
assemble into supramolecular fibrils with a precise axial 
stagger between neighbors of 67 nm, known as the D-period 
(6) (Fig. 1, C and D). Because each molecule spans a nonin¬ 
tegral number of periods (4.46 D), each period contains an 
overlap and a gap region that cause a cross-striated appear¬ 
ance in electron micrographs (Fig. 1 E). The D-banded axial 
periodicity is essential to ensure a high tensile strength (7,8), 
to regulate cellular recognition and binding of extracellular 
matrix molecules (9), and to template bone mineralization 
(10). Accordingly, mutations in type I collagen that cause 
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defects in assembly are associated with severe connective 
tissue disorders (11). 

In vitro assembly studies have provided key insights into 
the mechanism by which D-staggered collagen fibrils form. 
It was shown early on that purified collagen can form fibrils 
with the same axial periodicity as native fibrils without any 
accessory extracellular matrix components (12). Similar to 
native fibrils, the fibrils are tens of micrometers in length 
and -100 nm in diameter (13). Thus, all the required struc¬ 
tural information is encoded in the primary sequence. The 
kinetics of fibril assembly has been studied extensively, usu¬ 
ally by measuring the increase in solution turbidity that 
accompanies fibril formation. Turbidimetric growth curves 
of collagen exhibit an initial lag phase, followed by a growth 
phase during which the turbidity increases in a sigmoidal 
fashion (13-16). This behavior is characteristic of a coop¬ 
erative nucleation and growth mechanism. However, uncer¬ 
tainty remains about the exact nature of the nucleus and the 
degree of cooperativity (17—21) and about the mechanism of 
growth (22-25). There is some evidence that five-stranded 
microfibrils are already formed during the turbidimetric 
lag phase, which subsequently grow linearly and laterally 
during the growth phase (16,18). Several studies indicate 
that collagen monomers undergo conformational changes 
during the lag phase (23,26-28). Indeed, tropocollagen has 
to convert from a fluctuating, coiled conformation in solu¬ 
tion (Fig. 1 B) to a straightened conformation inside fibrils 

(Fig. 1 C). 

Collagen assembly is driven by noncovalent interactions, 
which are reversible upon changes in temperature or solution 
conditions. Assembly is promoted by warming and reversed 
by cooling, and is thus an endothermic process, similar to the 
assembly of other biopolymers such as actin and microtu¬ 
bules (29,30) and opposite to typical synthetic supramolecu¬ 
lar materials (31). Assembly is thus driven by an increase in 
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FIGURE 1 (color online) Schematic representa¬ 

tion of the hierarchical self-assembly of type I 
collagen. (A) Collagen molecules consist of 
three polypeptide chains twisted into a rod-like 
triple helix (not shown to scale; in reality 
the length is 300 and the diameter 1.5 nm). (. B ) 
The monomers are flexible polymers with a 
persistence length of ~10 nm. Shown is a 
typical conformation obtained by Monte Carlo 
simulations (courtesy of H. Amuasi). (C) Pen- 
tameric aggregate of collagen monomers in 
a D-staggered parallel arrangement. Note the 
straightened conformation of the collagen 
molecules. Each tropocollagen spans 4.46 D 
periods, giving rise to an overlap and a gap 
region within each period. ( D ) Fibrils are 
constructed from many pentamers. ( E) Transmission electron micrograph (2 /im wide) of negatively stained collagen fibrils assembled at 37°C in 
a physiological buffer clearly shows D-banding with an axial periodicity of 67 nm. 



entropy, resulting from the release of water molecules from 
the surface of tropocollagen helices (32). Opinions differ as 
to the relative contribution of different classes of noncovalent 
interactions to the thermodynamic driving force for assem¬ 
bly. Hydrophobic interactions are likely to play an important 
role. Atomistic simulations showed that water molecules are 
released from hydrophobic amino acid residues when the 
molecules pack into a fibril (33). Moreover, analysis of the 
amino acid sequence of collagen showed that hydrophobic 
residues occur in clusters and that D-staggered packing opti¬ 
mizes the alignment of hydrophobic side chains between 
apposing tropocollagens (34,35). Yet, there is also an alter¬ 
native interpretation of the endothermic nature of collagen 
assembly, namely that water is released from polar residues. 
According to this interpretation, the dominant driving force is 
provided by hydration forces associated with water-mediated 
hydrogen bonding between polar residues (36-38). This 
interpretation is supported by observations that sugars and 
polyols strongly reduce the attraction between collagen 
helices and suppress fibrillogenesis (19,39,40). These com¬ 
pounds are expected to disrupt water-mediated as well as 
direct hydrogen bonding between polar residues. Recent sim¬ 
ulations indeed indicate that clusters of hydrogen-bonded 
water exist in collagen fibrils, but the relative importance 
of hydration forces and hydrophobic forces could not be 
determined (33). The dependence of fibril assembly on solu¬ 
tion pH and ionic strength suggests that electrostatic inter¬ 
actions also reinforce the stability and/or axial specificity 
of assembly (20,41,42). The D-staggered packing arrange¬ 
ment appears to optimize alignment between complementary 
charges on opposing helices (34). Some studies also impli¬ 
cate specific ions such as divalent phosphate in providing 
salt bridges and improving axial order (15,42,43). Finally, 
there is evidence that the telopeptides further enhance the 
specificity of molecular recognition (44). When the telopep¬ 
tides are removed by pepsin or other proteases, nucleation is 
slowed down and the activation energy is increased (45—47). 
Computational models show that telopeptides dock to spe¬ 
cific helix recognition sites and induce reciprocal conforma¬ 


tional changes (48,49). It remains unknown whether this 
interaction is only catalytic in the initial stages of assembly 
or contributes also to the driving force. 

Despite several decades of intensive research on the 
process of collagen assembly, the molecular basis of the 
interactions responsible for collagen fibril stability (and 
mechanical strength) remains poorly understood. Here, we 
take an alternative approach to elucidate the role of nonco¬ 
valent interactions mediated by the collagen triple helix in 
fibrillogenesis. Rather than studying the process of assem¬ 
bly, we focus on the kinetics of cold-induced fibril disas¬ 
sembly. As a model system, we use collagen whose 
telopeptides are removed by pepsin. This is necessary to 
prevent spontaneous formation of intermolecular covalent 
cross-links (17,46,50-52). We used several complementary 
methods to characterize the kinetics of cold-induced disas¬ 
sembly of collagen fibrils formed at body temperature 
(37°C): we used turbimetry to probe the fraction of collagen 
in fibril form, atomic force microscopy (AFM) to probe the 
fibril morphology, and confocal microscopy and rheometry 
to measure the network connectivity. We demonstrate that 
collagen fibrils disassemble by subunit release from their 
sides as well as their ends, with complex kinetics involving 
an initial fast release followed by a slow release. Despite the 
absence of telopeptides, the fibrils are gradually stabilized 
over time. This dynamic stabilization may reflect structural 
plasticity of the collagen fibrils arising from their complex 
structure. In addition, we propose that the polymeric nature 
of collagen monomers may lead to slow kinetics of subunit 
desorption from the fibril surface. Dynamic stabilization 
of fibrils may be relevant in the initial stages of collagen 
assembly in vivo (2). 

MATERIALS AND METHODS 
Preparation of collagen gels 

Bovine dermal collagen I (NUTRAGEN) was purchased as a 6.0 mg/ml 
solution in 0.01 M HC1 from Advanced Biomatrix (San Diego, CA). This 
collagen is isolated from bovine hides by enzymatic digestion with pepsin. 
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which removes portions from the N- and C-terminal telopeptides containing 
the lysine groups that are responsible for intermolecular cross-linking in vivo 
(53). Sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS- 
PAGE) on 4% gels was used to assess the purity and degree of intramolecular 
cross-linking. Gels were stained with Instant Blue Coomassie staining solu¬ 
tion (Expedeon, Harston, UK). Samples with a final volume of 500 fiL and 
collagen concentration of 1 mg/ml were prepared by diluting the stock solu¬ 
tion with a buffer mixture, to give a final buffer composition of 10 mM phos¬ 
phate and 160 mM NaCl at pH 7.2 (I = 0.17). This is a near-physiological 
buffer that promotes a native D-banding pattern. In the case of turbidity 
experiments, the buffer mixture was degassed in a vacuum before mixing 
with collagen, to prevent the formation of air bubbles. All preparation steps 
involving collagen were performed at 4°C. Unless otherwise noted, chemi¬ 
cals were purchased from Sigma Aldrich (St. Louis, MO). The extent of 
collagen polymerization at 37°C was tested by precipitating the fibrils by 
centrifugation for 10 min at 16,100 x g in a 5415R tabletop centrifuge 
(Eppendorf, Hamburg, Germany) at 37°C. The collagen concentration in 
the supernatant was measured using the sirius red assay (54). The supernatant 
(100 fi\) was incubated for 30 min with 1 ml of 1 mg/ml Direct Red 80 in 1% 
picric acid. The solution was centrifuged for 10 min at 16,100 x g, the super¬ 
natant was removed, and the pellet was washed (to remove excess dye) with 
750 fi\ of 0.01M HC1. The pellet was finally dissolved in 250 fi\ 0.5M NaOH. 
After 5 min, the solution absorbance at 550 nm was measured using a 
Nanodrop 2000 spectrophotometer (Thermo Scientific, Waltham, MA) and 
calibrated against stock collagen in the same buffer. 


Turbidity measurements 

The kinetics of collagen polymerization at 37°C and subsequent depolymer¬ 
ization induced by cooling were probed by measuring the change in solution 
absorbance at a wavelength of 370 nm (A 370 ) with a Lambda 35 dual-beam 
spectrophotometer (PerkinElmer, Waltham, MA) equipped with a PTP 
6+6 Peltier temperature controller. Solutions of 1 mg/ml collagen were 
pipetted into a quartz cuvette (Quartz Suprasil precision cells 115B-S, 
Hellma Analytics, Mullheim, Germany) equilibrated to 37°C and absor¬ 
bance readings were initiated within 15 s. The baseline value was taken as 
the first recorded turbidity value. The absorbance readings were converted 
into turbidity values, r, by using the relation r = (A 370 In (10))/L, where L 
is the optical path length (1.0 cm). It was previously shown by comparison 
with fibril precipitation experiments that the turbidity of collagen solutions 
is directly proportional to the amount of mass present in fibril form 
(13,15). Thus, the turbidity provides a reliable, and model-independent, mea¬ 
sure of the relative degree of polymerization as a function of temperature. In 
some experiments we measured the wavelength dependence of the turbidity 
with 2 min time intervals; using light scattering theory modeling the fibers as 
rigid rods (55,56) it is possible to infer the fibril diameter and mass-length 
ratio (see the Supporting Material Note SI and Fig. SI). Collagen fibrils 
were always first allowed to assemble for 2 h at 37°C, and the absorbance 
was recorded with 30 s intervals. Subsequently, fibril disassembly was 
induced by lowering the temperature to target values between 4 and 37°C. 
We used two different disassembly protocols: temperature step experiments 
and temperature ramp experiments. The actual time needed for the sample to 
reach a given target temperature in step experiments varied from 1 to 4 min, 
depending on the target temperature, as measured with a thermocouple. In 
temperature ramp experiments, the temperature was lowered at a constant 
rate between 0.2 and 129°C/h. The actual time dependence of the tempera¬ 
ture was measured with a thermocouple and used in the analysis. The solution 
absorbance during disassembly was measured for 2 h at 30 s intervals in step 
experiments and at 15 min intervals during ramp experiments. 


AFM and transmission electron microscopy (TEM) 

To characterize fibril morphology and diameter, we imaged fibrils deposited 
on Formvar-coated glass coverslips using AFM (Dimension 3100 AFM, 


Veeco digital instruments, Plainview, NY). We were unable to image the 
fibrils in buffer, because the fibrils are very soft and easily pollute the 
AFM tip or dislodge from the surface. For this reason, we imaged the fibrils 
in air. For comparison, we also performed TEM at 80kV (Tecnai G2, FEI 
company, Hillsboro, OR) on samples that were prepared identically, but 
on Formvar-coated copper electron microscopy (EM)-grids. The glass sub¬ 
strates for AFM were cleaned with 70% ethanol and dried with a flow of N 2 , 
and then briefly dipped into a Formvar solution (1% in 1,2-dichloroethane) 
and slowly extracted to ensure a homogeneous layer. The Formvar was 
allowed to dry in air. The substrate was placed gently for 10 s on top of a 
collagen gel (polymerized in a humified petri dish) at the appropriate tem¬ 
perature. The substrate was washed 3 times with phosphate buffered saline 
(PBS) of pH 7.2 and 3 times with milliQ, water to prevent the formation of 
salt crystals. All solutions were of the appropriate temperature and excess 
liquid was blotted off with filter paper after each washing step. The sub¬ 
strate was dried at 37°C for at least 2 h. The fibrils were imaged in tapping 
mode using a TESPA doped silicon cantilever (Bruker AXS, Camarillo, 
CA) with a nominal tip radius of 8 nm and spring constant of 42 N/m. 
Images were flattened using the Nanoscope software (6.14R1, Veeco digital 
instruments), but otherwise unaltered. The fibril width and height were 
analyzed in ImageJ software (version 1.46a, National Institutes of Health, 
Bethesda, MD), averaging over 20 fibrils per experimental condition. 

Confocal reflectance microscopy 

The three-dimensional architecture of collagen gels in their native, hydrated 
state was examined by confocal microscopy with an Eclipse Ti inverted 
microscope (Nikon instruments Europe, Amstelveen, The Netherlands). 
We used the reflectance mode, to avoid the use of fluorescent labels, which 
can influence collagen (de)polymerization (57-59), using 457 nm light 
from an Ar laser (Melles Griot, Albuquerque, NM). Collagen gels were 
polymerized for 2 h at 37°C in glass chambers made from a microscope 
slide and coverslip separated by Parafilm spacers and sealed with vacuum 
grease. Gels were imaged either directly after polymerization or after 
subsequent cooling to 4°C. Image stacks were obtained at least 10 fin 1 
away from the glass surface by scanning through the z direction in steps 
of 0.5 fim over a range of 20 fim with a piezo-driven lOOx oil immersion 
objective (Nikon). To optimize the signal/noise ratio, we used 2x line 
average. Maximum intensity z projections were made with ImageJ. 

Cone-plate rheology 

The linear viscoelastic properties of the collagen gels were measured using 
a stress-controlled MCR501 rheometer (Anton Paar, Graz, Austria) with a 
CP40-2 acrylic measuring cone with 40 mm diameter and 2° angle and a 
steel bottom plate heated to 37°C with a PTD200 Peltier system. The 
collagen solution was applied to the prewarmed rheometer, and allowed 
to polymerize in situ for 2 h before starting oscillatory shear measurements. 
The sample was in a closed hood to maintain a constant temperature and 
humidity. The frequency-dependent shear modulus G*(cu) = (t(oj)/ y(w) of 
the collagen gels was measured by applying an oscillatory strain y(w) 
decreasing logarithmically from 34 to 0.0628 rad/s, and recording the stress 
response, a(oj). The strain amplitude was 0.5%, which is much smaller than 
the critical strain (~5%) where nonlinearity sets in. G*(w) is a complex 
modulus with an in-phase (elastic) component, G'(gj), and an out-of-phase 
(viscous) component, 

RESULTS 

Kinetics of collagen assembly and disassembly 

Collagen self-assembly is entropy-driven and depends on 
noncovalent interactions (20,32). Fibrils can be assembled 
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by heating a collagen solution from 4 to 37°C and disas¬ 
sembled by subsequent cooling. We use turbidimetry to 
analyze the kinetics of the disassembly process in response 
to a sudden stepwise lowering of the temperature, because 
the solution turbidity is a direct measure of the amount of 
mass present in fibril form (13,15). Fig. 2 a shows three 
typical assembly-disassembly experiments, with identical 
assembly conditions (1 mg/ml collagen in PBS, 2 h at 
37°C) but differing in disassembly temperature (r dis = 4, 
17, 27°C). The change in turbidity during the assembly 
phase is highly reproducible, showing three distinct phases 
characteristic of a nucleation-and-growth process (14). 
There is an initial lag phase of -10 min with no turbidity 
change. Prior studies indicate that short, thin fibrils are 
already present at the end of the turbidimetric growth phase 
with -100 collagen molecules per unit length, but the struc¬ 
tures are too small during the lag phase to give a measurable 
turbidity (16,22,60). This is confirmed by our own wave- 
length-dependent turbidity measurements (Note SI and 
Fig. S2 a). The lag phase is followed by a growth phase 
with a sigmoidal increase in turbidity, during which the 
collagen fibrils grow in length and width. Wavelength- 
dependent turbidity data show that the mass-length ratio 
of the fibrils also increases sigmoidally during the growth 
phase (Fig. S2 a). Finally, all three turbidity curves reach 
a plateau phase where the turbidity reaches a constant level 
(2.57 ± 0.06), which varies by <4% among experiments. It 
actually takes 6 h for the turbidity to reach a constant value, 
and at 2 h, the turbidity has reached -93% of this maximum 
(consistent with precipitation experiments, showing that 
collagen assembly after 2 h is 87% complete). Nevertheless, 
for practical reasons we always assemble gels for only 2 h, 
and denote the fibrils at this point as mature, with a turbidity 
value T 37 . Wavelength-dependent turbidity data indicate that 
the mature fibers have an average diameter of 61 nm 
(Fig. S2 b) with 1070 monomers per cross section (Note 
SI), consistent with prior reports for rat tail collagen (61). 

After 2 h, we reduce the temperature stepwise, which 
causes an immediate drop in the turbidity (Fig. 2 A). The 
slope of the turbidity response curve is relatively shallow 
in the first few minutes, because of the finite duration of 
the temperature step (1-4 minutes). The turbidity value 


r d i S reached after 2 h of disassembly is strongly dependent 
on the disassembly temperature: the lower the temperature, 
the lower r dis , indicating that the fibrils release more mass. 
The kinetics of disassembly clearly does not follow a simple 
single exponential form. Empirically, the decay curves are 
well fit by a sum of two exponentials, corresponding to a 
fast decay with a characteristic decay time of -9 min, inde¬ 
pendent of disassembly temperature, followed by a slow 
decay with a decay time that increases from 200 min at 
4°C to 2000 min at 27°C (Fig. S3). At 32°C, the turbidity 
is essentially constant at long times. As shown in Fig. S4, 
the turbidity during a disassembly experiment at 22°C is 
actually still not constant after 4 days. 

To quantify the extent of disassembly, we take the ratio 
between the turbidity measured after disassembly and that 
measured for the mature fibrils, T dis /T 37 which is a measure 
of the degree of fibril stability. For practical reasons, we mea¬ 
sure this ratio after 2 h of disassembly. Although the disas¬ 
sembly process is <80% complete at this time (Fig. S4), 
the rate of disassembly is already much smaller than directly 
after the temperature step, making this an acceptable choice. 
As shown in Fig. 2 B, T dis /T 37 decreases nearly linearly when 
the disassembly temperature is lowered from 32 to 7°C. The 
turbidity does not return to zero upon cooling, even to 4°C, 
which indicates that the collagen fibrils do not completely 
disassemble. To test whether collagen fibrils indeed remain 
in solution, we cycle samples cooled down to 22° C for 2 h 
back to 37°C. As shown in Fig. S5, disassembly is reversible, 
because heating to 37°C brings the turbidity back to its 
original value. Strikingly, there is no lag phase during the 
reassembly process, which strongly indicates that collagen 
fibrils indeed remain and can act as templates for fibril growth 
upon reheating. The wavelength dependence of the turbidity 
indicates that fibrils disassembled for 2 h at 4°C still contain 
229 monomers per cross section, which is -20% of the mass/ 
length ratio of mature fibrils (Fig. S2 a). The loss in mass 
from the sides (/u dis //x 37 ) is comparable to the overall loss 
in fibril mass (T dis /r 37 ) at all disassembly temperatures 
(Fig. S 6 ), suggesting that fibrils predominantly shed mono¬ 
mers from their sides. 

The apparent stability of the fibrils at 4°C is surprising, 
because collagen solutions kept at 4°C and never exposed 



time (hours) disassembly temperature(°C) 


FIGURE 2 Turbidity measurement of the assem¬ 
bly and disassembly kinetics of a 1 mg/ml collagen 
solution, (a) Turbidity response for three collagen 
samples, each assembled for 2 h at 37°C and disas¬ 
sembled by cooling to different temperatures, Tdis- 
( b ) Temperature dependence of TdijT?>i (note the 
inverted temperature axis). 
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to 37°C do not develop any measurable turbidity (Fig. S7). 
This history dependence suggests that fibrils formed at 37°C 
and cooled to 4°C are kinetically stabilized. It should be 
noted that, although a collagen solution at 4°C does not pro¬ 
duce aggregates detectable by the turbidity method, some 
aggregation does take place at 4°C. When a neutralized 
collagen solution is heated to 37°C after 17 h incubation 
at 4°C, the resulting turbidity curve has no lag phase (inset 
of Fig. S7). Apparently, there is a nonzero driving force for 
collagen to self-assemble at 4°C into structures that can act 
as nuclei for fibril growth. The rate of assembly after 17 h 
incubation at 4°C is slower than the rate of assembly of a 
freshly neutralized collagen solution (cf. Fig. 2 A), suggest¬ 
ing that the oligomeric structures formed at 4°C are different 
from the aggregates formed in the nucleation phase at 37°C. 
We conclude that there is a clear history dependence 
of the turbidity measured at 4°C, suggestive of kinetic 
stabilization. 

To test whether the stability of collagen fibrils at low tem¬ 
perature is truly kinetic in origin, we examined whether 
collagen gels could be redissolved in 0.01 M HC1. Under 
these conditions, there is no driving force for collagen fibril- 
logenesis (62). We polymerize collagen at 1 mg/ml for 2 h 
and then dialyze against 0.01M HC1 at 4°C. To test whether 
fibril disassembly is complete, we perform an assembly 
experiment at 37°C with the dialyzed solution. If the fibrils 
were fully dissolved, the assembly process should have a lag 
phase, similar to the first assembly step in Fig. S5. However, 
if parts of fibrils or other nuclei remain because they are 
being held together by covalent cross-links, there should 
be no lag phase. For the assembly experiment we use a con¬ 
centration around 0.1 mg/ml, where the lag phase is longer 
and easier to recognize than at 1 mg/ml (13). We indeed find 
that the lag phase returns (Fig. S8), which implies that the 
collagen fibrils are in principle fully soluble. This indicates 
that there is some kinetic mechanism that stabilizes the 
fibrils in a near-physiological buffer at neutral pH. More¬ 
over, this stabilization is apparently an intrinsic property 
of the collagen triple helix, because we use pepsin-treated 
collagen that lacks the telopeptides. 

Rate dependence of fibril disassembly 

The experiments described previously show the response of 
the system to an abrupt change of temperature. This is a 
nonequilibrium experiment, because the temperature is 
changed much faster than the equilibration time of fibril 
disassembly. To obtain more insight into the fibrils’ equili¬ 
bration behavior, we perform a series of experiments where 
the temperature is lowered from 37 to 22°C at a constant 
rate. Because the temperature jump experiments show a 
fast disassembly process with a decay time of 9 min and a 
slow process with a time constant of several days, we choose 
a range of cooling rates from 0.2°C/h (which takes over 
3 days) to 128.6°C/h (which takes 7 min). Fig. 3 shows a 



time (hours) 


FIGURE 3 Normalized turbidity response (t/t 37 ) of a 1 mg/ml collagen 
solution during a temperature ramp experiment. There are three stages, as 
indicated by the dotted line showing the temperature: 1), an assembly stage 
for 2 h at 37°C, 2), a temperature ramp stage in which the temperature is 
lowered to 22°C at a constant rate, and 3), a tail stage, during which the tem¬ 
perature is maintained at 22°C. 


typical turbidity response curve, which encompasses three 
stages, as indicated by the dotted line showing the tempera¬ 
ture: 1), assembly for 2 h at 37°C, 2), temperature ramp 
down to 22°C with a constant rate, and 3), a tail stage at a 
constant temperature of 22°C. In the assembly stage, fibrils 
are formed, producing a sigmoidal increase in turbidity. 
Next, the temperature is ramped down toward 22°C at a 
rate of 5°C/h. This does not lead to an immediate change 
in t/t 37 likely because the fibrils are not yet fully assembled 
after 2 h, and continue to grow as long as the temperature is 
sufficiently close to 37°C. After a short delay, though, t/t 37 
drops at a constant rate until the target temperature is 
reached. Strikingly, t/t 37 continues to decrease at 22°C, at 
a rate that gradually slows with time. We call this stage 
the tail stage. 

Remarkably, the fibrils’ response during the temperature 
ramp stage is independent of the cooling rate within exper¬ 
imental error (Fig. 4). Note that, because the sampling rate is 
limited to one measurement every 15 min, the faster ramps 
have only a few data points. Nevertheless, these data clearly 
show that t/t 37 is a function of the temperature alone, and 
not of the rate of change of the temperature. This is true 
for a broad range of ramp rates, with the fastest rate being 
over 600 times faster than the slowest one. This result sug¬ 
gests that a certain fraction of collagen monomers on the 
surface of the fibrils is quickly exchangeable. However, 
this is clearly not true for all monomers, because after the 
ramp stage, the fibrils continue to disassemble at a slow 
rate. As shown in Fig. 5 A, the extent of disassembly at 
22°C depends on the cooling rate of the preceding ramp 
stage. The higher the cooling rate during the ramp, the larger 
the decrease of t/t 37 . The disassembly curves have a com¬ 
plex functional form, which we did not attempt to fit. The 
initial rate of disassembly, measured as the initial slope of 
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FIGURE 4 {color online) Normalized turbidity response (T djs /T 37 ) of five 
different collagen networks during the temperature ramp stage of tempera¬ 
ture ramp experiments. In each case, the temperature is lowered from 37°C 
to 22°C, but at different rates, as indicated in the legend. 

t/t 37 , increases approximately linearly with the cooling rate 
(Fig. 5 B). This thermal memory effect is strongly sugges¬ 
tive of a slow, time-dependent stabilization of the collagen 
fibrils: when the fibrils are cooled at a slower rate, there is 
more time for fibril stabilization. To test this hypothesis, 
we assemble collagen fibrils at 37°C for various lengths of 
time and measure the turbidity change in response to abrupt 
cooling to 22°C. In accordance with our hypothesis, we 
observe that the apparent stabilization (expressed in terms 
of T dis /T 37 ) increases from T dis /T 37 = 0.24 to 0.75 when the 
incubation time at 37°C is increased from 1 to 18 h 
(Fig. S9). 

Disassembly causes fibril thinning but leaves the 
fibril structure intact 

Turbidimetry provides an ensemble measure of the degree 
of fibril disassembly, but gives no information on the 
changes in fibril morphology or network structure upon 
cooling. It is conceivable that disassembly might change 
the packing structure of the fibrils. To directly visualize 
effects of disassembly on fibril morphology, we use AFM 
to image collagen fibrils directly after assembly (2 h at 
37°C) and after disassembly at either 22 or 4°C for 2 h. 


As shown in Fig. 6, mature collagen fibrils are long and 
straight and have an average apparent width of 221 ± 
28 nm and height of 30 ± 9.6 nm. This flattening was 
also observed elsewhere (63), and is most likely a conse¬ 
quence of adsorption to the substrate and (partial) loss of 
hydration water upon drying in air (64). The measured fibril 
width and height are therefore not directly representative of 
the diameter of fibrils in their native, hydrated state. How¬ 
ever, we can still use AFM to measure relative differences 
among different temperature conditions. In most of the 
fibrils the distinctive 67 nm D-periodic banding is visible, 
as exemplified by the axial height profile shown underneath 
the image. 

Disassembly at 22 or 4°C for 2 h clearly changes the fibril 
dimensions. The average apparent width decreases to 175 ± 
25 nm at 22°C and 146 ± 34 nm at 4°C, whereas the 
average height decreases to 22 nm in both cases (± 9 nm 
at 22°C and ± 4 nm at 4°C). These changes indicate that 
the fibrils lose mass from their sides, consistent with the 
reduction of the fibril mass-length ratio measured by wave- 
length-dependent turbidimetry (Fig. S6). However, the 
major morphological features are not visibly affected by 
disassembly. The fibrils still appear long and straight, and 
retain their D-banding, indicating that the native intermolec- 
ular stagger is still intact. The fibrils still appear as flattened 
cylinders, structurally the same as mature fibrils, which sug¬ 
gests that no large groups of monomers break off during the 
disassembly process. We conclude that disassembly occurs 
all along the fibril by dissociation of monomers or small 
aggregates from the outer surface of the fibril, with the re¬ 
maining fibril core retaining its structure and organization. 

A striking feature of the AFM data is the material in the 
background. For the mature fibrils there is a small amount 
of material, consisting of amorphous aggregates. For the 
partially disassembled fibrils, the background material is 
more abundant, especially for the fibrils disassembled at 
4°C. For fibrils disassembled at 4°C the background mate¬ 
rial looks fibrillar. However, it is uncertain whether the 
background material is also present in solution. We cannot 
rule out that it may result from the AFM sample preparation. 
The Formvar-coated glass substrate has a high affinity for 
collagen, and may therefore cause monomers or oligomers 
to aggregate on the surface (65). Similar aggregates appear 



time after completion of temperature ramp (hours) 



FIGURE 5 (Left) Normalized turbidity response 
(Tdi S /r 37 ) of five different collagen networks during 
the tail stage of temperature ramp experiments. For 
a larger preceding ramp rate, there is more disas¬ 
sembly during the tail stage. (Right) The initial 
slope of each tail response curve is a linear function 
of the preceding ramp rate. 
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FIGURE 6 Atomic force micrographs of collagen fibrils assembled at 1 mg/ml, (a) directly after assembly or after subsequent 2 h disassembly at (b) 22°C 
or (c) 4°C. Representative height profiles are measured along a fibril in each image (white lines in the insets), showing that the 67 nm D-periodicity is present 
under all three conditions. Insets (0.78 X 0.78 /ini 2 ) show an enlarged view of the measured fibrils. 


in TEM images of fibrils on Formvar-coated electron micro¬ 
scopy (EM) grids (Fig. S10). 

Disassembly minimally affects network 
connectivity 

The AFM images demonstrate that the turbidity decrease 
upon cooling is, at least partly, caused by monomer loss 
from the sides. However, they do not reveal any effects on 
fibril length because the fibrils are many micrometers long 
(22). To test whether fibrils also lose mass from their 
ends, we probe the change in network connectivity at 4°C 
with two complementary methods. We use confocal reflec¬ 
tance microscopy (CRM) to visualize the networks in their 
native, hydrated state and reveal their connectivity, homoge¬ 
neity, and mesh size (57). Moreover, we use cone-plate 
rheology to measure the elastic modulus of the gels, which 
is sensitive to network connectivity (66). 

Fig. 7 shows two maximum intensity projections of 
20 /xm-thick sections of collagen gels, one taken of a gel 
formed at 37°C for 2 h (A) and the other taken after subse¬ 
quent disassembly at 4°C for 2 h (B). At first sight, both 


images look indistinguishable: the networks are homoge¬ 
neous and isotropic and have a similar fibril density. In 
contrast to the AFM images, no change in fibril diameter 
can be observed. Indeed, the fibrils have diameters below 
the diffraction limit, so a diameter reduction is only ex¬ 
pected to show up in a CRM image when it is so extreme 
that the fibrils do not backscatter any light. In the mature 
gel at 37°C, there are no visible fibril ends and also time 
lapse imaging shows that the fibrils are motionless. Thus, 
we are unable to determine the lengths of the fibrils. In 
the gel that has been disassembled at 4°C we can still 
not recognize fibril ends in still images, but time lapse 
imaging reveals some dangling fibril ends as recognized 
by their transverse fluctuations. This effect is illustrated in 
Fig. 7 C, which shows five images with an undulating fibril 
end in the middle, recorded 4.75 s apart. The position and 
shape of the fibril in the first frame are marked with five 
white dots in each frame, which clearly reveals motion of 
the fibril end. The maximum intensity projection (panel 
all) shows that the fibril end undulates, whereas the other 
fibrils are motionless. Apparently, there is some endwise 
shrinkage of filaments. However, the loss in fibril mass 



FIGURE 7 Maximum intensity projections of 
z-stacks (41 z planes spaced by 0.5 /xm) of a 
1 mg/ml collagen gel, (a) directly after 2 h of assem¬ 
bly at 37°C, and ( b ) after subsequent disassembly 
for 2 h at 4°C. The bright spot in the center of the 
images is an artifact of CRM caused by reflection, 
(c) Five frames of a time lapse movie (taken 4.75 
seconds apart) of a 1 mg/ml collagen gel after partial 
disassembly at 4°C for 2 h. The dotted line indicates 
a fibril with a dangling end. The last image (all) 
shows a maximum intensity projection that clearly 
reveals fluctuations of the dangling fibril end. 
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(85% at 4°C) can clearly not be explained by endwise 
shrinkage alone, because then we should see only 15% of 
the original fibril density in the confocal images of disas¬ 
sembled gels. Indeed, the turbidity data show that the overall 
loss in fibril mass is close to the loss in mass-length ratio 
(Fig. S6): at 4°C, the overall mass loss relative to mature 
fibrils is 85%, whereas the loss from the sides is 79%. 

To verify that disassembly indeed does not disrupt 
network connectivity, we also performed rheology on gels 
before and after disassembly. Fig. 8 shows that after assem¬ 
bly for 2 h at 37°C, the storage modulus, G', is 30 Pa 
at 1 rad/s, consistent with prior reports (57,66). The gels 
behave as weak elastic solids, with a slight frequency depen¬ 
dence of G' and a loss tangent G"/G' of 0.25 (inset). After 
subsequent disassembly for 2 h at 4°C, the G' is reduced 
~70-fold at an angular frequency of 10 rad/s. However, 
the sample is still a solid, with an elastic modulus that is 
only weakly dependent on frequency and a loss tangent 
that is still below 1 (0.08). Despite losing 85% of its mass, 
the fibril network remains space-spanning and elastic. The 
reduction of G' may be caused by the increase in elastically 
inactive dangling ends observed by confocal microscopy. 
Prior rheology measurements of collagen networks have 
shown that G' depends strongly on protein concentration, 
usually scaling as a power law in concentration with an 
exponent between 2 and 3 (61,66,67). In addition, the reduc¬ 
tion of G' may reflect the reduced fibril rigidity caused 
by the substantial reduction in fibril diameter. There 
is not yet a quantitative theoretical model that predicts 
the modulus of collagen networks. There is evidence that 
collagen networks deform in a highly nonaffine (nonuni¬ 
form) manner (68). Nonaffinity is known to make networks 
much softer than they would be in the affine limit, but there 
is no quantitative (analytical) prediction of the extent of 



FIGURE 8 Network storage ( elastic ) modulus as a function of angular 
frequency, for a 1 mg/ml collagen gel subjected to a temperature cycle: 
fibrils were first assembled for 2 h at 37°C, then cooled down for 2 h at 
4°C, and finally reassembled for 2 h at 37°C. The inset shows the corre¬ 
sponding loss tangents. 


softening (69). When we reassemble the gel for 2 h at 
37°C, G' and G" both return close to their original values. 
This is in accordance with the turbidity assay showing a 
return of the original solution turbidity (Fig. S2). 

DISCUSSION 

Prior studies of collagen reversibility focused on covalent 
stabilization by the telopeptides, which is crucial for the 
tensile strength of collagen fibrils in vivo. Assembly of 
neutral-salt soluble or acid-solublized collagen, both of 
which possess intact telopeptides, is only partly reversible 
(52,70). Already during polymerization, allysine residues 
in the telopeptide regions condense with lysine or hydroxy- 
lysine residues in helix domains of adjacent molecules in 
the fibril, creating intermolecular cross-links. Upon aging 
of collagen gels at 37°C, assembly becomes progressively 
more irreversible (50,71). When cross-linking is prevented 
by using lathyritic collagen from animals fed with an 
inhibitor of lysyl oxidase, assembly is more reversible 
(72). Yet, there are a few prior turbidimetric reports hinting 
at a time-dependent increase in stability of fibrils formed 
from pepsin-collagen (73,74). The buffer composition was 
different from ours (pH 6.8, 300 mM phosphate), but the 
observations were comparable. In response to a sudden 
temperature drop, the turbidity decayed in a biphasic 
manner and gels aged at 37°C became progressively more 
insoluble. However, no imaging data were shown, so it is 
unclear what was the packing arrangement of monomers 
(which may not be D-staggered at 300 mM phosphate) or 
the relative contributions of endwise and sidewise disas¬ 
sembly. Similar observations were reported for lathyritic 
collagen, which also showed thermal memory effects 
despite the absence of cross-linking (72). 

What is the mechanism of the time-dependent stabiliza¬ 
tion of fibrils? One viable mechanism is structural plasticity. 
This phenomenon was proposed to cause dynamic stabiliza¬ 
tion of two other supramolecular protein polymers, actin 
filaments and microtubules (75,76). Actin filaments, like 
collagen fibrils, display a biphasic disassembly response 
and a progressive stabilization upon aging. Based upon 
EM and x-ray scattering data, it was proposed that filaments 
can gradually change from a more disordered and unstable 
structural state into a more ordered and stable structural 
state. Recently, however, it was proposed that the reported 
stabilization of actin was due to photo-induced cross-linking 
induced by the intense illumination employed in the fluores¬ 
cence microscopy assays used for monitoring disassembly 
(77). In our study, we robustly find time-dependent stabiliza¬ 
tion using techniques which do not employ illumination 
(EM, AFM, rheology). We propose that collagen fibrils 
may be relatively disordered directly after assembly, but 
acquire more order over time, by internal rearrangements. 
It is a formidable challenge, however, to directly demon¬ 
strate structural plasticity for collagen fibrils, given their 
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enormously complex and semicrystalline three-dimensional 
architecture. While actin filaments and microtubules 
behave as linear polymers whose assembly dynamics are 
restricted to the ends, collagen fibrils have hundreds of 
monomers per cross section packed in an array whose lateral 
packing order remains poorly understood. Our EM and AFM 
images indicate that axial order is established immediately, 
probably because of the high degree of specificity of the 
noncovalent interactions that drive assembly. However, we 
have no information on the lateral packing order. X-ray 
diffraction measurements of native collagen fibrils show 
varying degrees of lateral disorder (6,78). Fibrils assembled 
in vitro were shown to have an even less well-ordered lateral 
packing (12). Synchrotron x-ray scattering measurements 
or EM microscopy of fibril cross sections may shed more 
light on time-dependent changes in lateral packing (79). 
Alternatively, antibodies might be developed that recognize 
specific structural states of collagen fibrils in analogy to actin 
and microtubules (80,81). A complementary approach may 
be computational modeling of collagen fibril assembly using 
coarse-graining approaches (82), possibly combined with 
atomistic simulations (5,33). 

There are additional complexities in collagen assembly 
that may contribute to kinetic stabilization of the fibrils. 
Unlike actin and tubulin monomers, which are compact 
globular subunits, collagen monomers are long polymers 
with substantial flexibility (4,5). Both monomer additions 
to the surface of growing fibrils and monomer removal 
from disassembling fibrils require large conformational 
changes of the monomers, which may pose large kinetic 
barriers. There is indeed some evidence from hydrody¬ 
namic assays and vibrational spectroscopy of intramolecular 
conformational changes during fibril nucleation and growth 
(23,26-28). A tractable way to model these processes may 
be by using theories of polymer crystallization and desorp¬ 
tion (83-85). Desorption of semiflexible polymers from a 
surface is known to be slow. Moreover, it is possible that 
collagen disassembly sets up a complex ecosystem of 
different supramolecular species that can exchange with 
the sides and ends of fibrils. For instance, oligomers may 
be released in addition to monomers and released monomers 
may potentially associate in solution. 

Dynamic stabilization of fibrils may be relevant in the 
initial stages of collagen assembly in vivo, particularly in 
growing tissues and during fibrosis and wound healing (2). 
It will be interesting to test whether other extracellular 
matrix compounds such as proteoglycans affect the kinetics 
of cold-induced disassembly of collagen (3). Our results 
may be of even more direct relevance in the context of 
biomedical applications. One important area is tissue optical 
clearing for biomedical optics and photomedicine (86). It is 
known that the turbidity of biological tissues can be reduced 
by immersion in chemical agents such as sugar alcohols and 
restored by immersion in physiological saline. It was pro¬ 
posed that these agents reversibly weaken noncovalent bind¬ 


ing forces. Our results may help to explain the basis of this 
reversibility. A second application area where our results are 
relevant is in tissue regeneration and drug delivery (87). It is 
crucial to understand the molecular basis of collagen stabil¬ 
ity to control the kinetics of in vivo remodeling, resorption, 
and controlled drug release. 

CONCLUSIONS 

We demonstrate that collagen fibrils disassemble by subunit 
release from their sides as well as their ends, with complex 
kinetics involving an initial fast release followed by a slow 
release. Despite the absence of telopeptides, the fibrils are 
gradually stabilized over time. On a timescale of several 
days, the fibrils are stable even at temperatures as low as 
4°C. Yet, when collagen is never heated and kept at 4°C, there 
is no formation of fibrils. This is clearly indicative of a kinetic 
stabilization. Thermal memory is also evident from temper¬ 
ature ramp experiments: when collagen solutions are cooled 
at a constant rate, the turbidity decreases in a rate-indepen- 
dent manner during cooling, but continues to decrease slowly 
at a rate that is correlated with the preceding cooling rate once 
the temperature has reached a constant value. Comple¬ 
mentary experiments show that the stability of collagen 
fibrils formed at 37°C depends on age. We can rule out any 
irreversible changes such as denaturation, because the fibrils 
are perfectly soluble in acid and fibrils can be reverted to their 
original state by reheating to 37°C. We propose that the com¬ 
plex, hierarchical structure of collagen fibrils is responsible 
for the kinetic stabilization. The dynamic stabilization may 
reflect structural plasticity of the collagen fibrils arising 
from their complex structure or slow kinetics of subunit 
desorption from the fibril surface related to the polymeric 
nature of the collagen subunits. Dynamic stabilization of 
fibrils may be relevant in the initial stages of collagen assem¬ 
bly during tissue morphogenesis, fibrosis, and wound healing. 

SUPPORTING MATERIAL 

Supporting data (10 figures with additional data and one Supporting Note 
explaining the determination of collagen fiber diameter and mass-length 
ratio by turbidimetry) are available at http://www.biophysj.org/biophysj/ 
supplemental/S0006-3495( 13)00619-X. 
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